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ABSTRACT 


The  solution  of  a  particular  electrostatic  potential  problem  in¬ 
volving  a  conducting  sphere  and  infinite  ground  plane  in  the  presence 
of  a  point  charge  is  discussed.  Use  is  made  of  the  conventional  im¬ 
aging  method  to  obtain  the  solution  in  terms  of  'Lambert  series .  It 
is  shown  how  some  of  the  resultant  series  may  be  combined  and  expressed 
as  contour  integrals  in  the  complex  plane,  and  how,  by  a  suitable  defor¬ 
mation  of  the  contour,  the  series  may  be  transformed  into  an  extremely 
rapidly  convergent  series .  The  contour  method  is  used  to  obtain  very 
simply  an  asymptotic  expansion  of  part  of  the  solution  for  a  limiting 
case  of  the  problem  which  was  obtained  very  laboriously  by  Schlomilch. 

1.  INTRODUCTION 

»-  ’  '  .  ,  ,  .  if  ,* 

The  problem  in  electrostatics  of  determining  charges  and  potentials 
on  conducting  spheres  and  ground  planes  in  the  presence  of  each  other 
and  in  various  combinations  has  been  of  interest  to  physicists  for  many 
years.  The  problems,  in  general,  are  most  easily  solved  by  the  method 
of  images,  to  which  there  are  many  standard  references  (ref  l).*  The 
solutions  for  the  charges  and  potentials  on  the  spheres  are  in  the  form 
of  sums  of  Lambert  series  (ref  2)  with  coefficients  an  of  the  form  \h/2, 

where  K  may  be  positive  or  negative. 

/ 

The  Lambert  series  had  been  expressed  in  several  forms;  however, 
none  were  good  for  computing  results  when  two  of  the  conducting  bodies 
are  very  close.  In  i860,  Schlomilch  (ref  j)  devised  a  method  which  was 
later  extended  by  Russell  (ref  by  which  the  very  close  case  can  be 
expressed  in  an  asymptotic  series  which  is  more  accurate  as  the  bodies 
move  closer  together. 

Even  with  the  contributions  of  Schlomilch  and  Russell,  there  was 
still  a  need  for  more  improvement  in  the  computational  methods.  The 
rapidly  convergent  forms  for  the  Lambert  series  have  individual  terms 
which  require  laborious  computation  in  themselves  and  for  this  reason 
are  no  more  desirable  for  computation  than  a  more  slowly  convergent 
Lambert  series  which  has  relatively  simple  individual  terms. 

In  many  problems,  two  of  the  Lambert  series  in  the  solution  com¬ 
bine  to  form  a  single  series  which  is  summed  from  -oo  to  +oo.  This  ser¬ 
ies  may  be  expressed  as  a  contour  integral  in.  the  complex  plant  (ref  5) • 

A  proper  adjustment  of  the  path  of  integration  yields  a  single  series 
which  is  so  rapidly  convergent  that  any  more  than  the  first  two  terms  are 
rarely  required. 

•  i  .  i  ;  c 

Below,  a  specific  potential  problem  involving  a  conducting  Sphere 
and  plane,  and  a  point  charge  is  handled  in  some  detail.  The  Standard 


*  References  are  listed  on  page  20. 
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imaging  methods  are  used  to  obtain  the  solution,  and  the  more  powerful 
methods  described  above  are  applied  to  reduce  the  solution  to  a  more 
usable  form. 

2.  THEORY 

In  the  theory  of  electrical  images  it  is  shown  that  a  conducting 
sphere  of  radius  a,  (figure  l)  in  the  presence  of  a  positive  point 
charge  q  located  at  distance  d  from  its  center  will  be  an  equipotential 
surface  if  a  negative  image  charge  of  magnitude  ^  is  placed  on  the 

axis  joining  the  charge  and  sphere  at  a  distance  — from  the  center  of 
the  sphere  (ref  l).  This  combination  causes  the  sphere's  surface  to  be 
at  zero  potential.  If  we  wish  the  sphere  to  be  at  any  other*  potential, 
we  may  place  an  additional  charge  at  the  center  so  that  the  potential 
of  the  sphere  Is  this  charge  divided  by  the  radius.  Similarly,  the  im¬ 
age  which  keeps  an  infinite  plane  an  equipotential  in  the  presence  of 
a  point  charge  is  a  charge  of  equal  magnitude  and  opposite  sign  placed 
an  equal  distance  below  the  plane. 


When  we  consider  the  problem  of  a  charged  sphere  and  a  plane  as 
shown  in  figure  2,  we  find  that  an  infinite  number  of  images  are  nec¬ 
essary  to  meet  the  conditions  that  the  plane  and  sphere  are  equipotential 
surfaces.  The  charge  g  at  the  center  of  the  sphere  images  in  the  plane 
as  -Q  to  keep  the  plane  at  equipotential.  Now,  however,  the  sphere  in 
the  presence  of  -Q,  is  not  at  equipotential  until  -Q  is  imaged  back  in  the 
sphere.  This  image  charge  makes  the  plane  nonequipotential  until  it 
images  back  in  the  plane,  and  so  the  process  continues.  The  result  is  an 
infinite  series  of  charge  images.  The  potential  of  the  sphere  is  deter¬ 
mined  by  the  charge  at  Its  center,  i.e.,  V  =  — ,  where  a  is  the  sphere 
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•  -Q-l 

•  -Q 


Figure  2 


radius,  and  the  total  charge  on  the  sphere  is  the  sum  of  the  center 
charge  and  all  the  image  charges  in  the  sphere.  The  plane,  of  course, 
stays  at  zero  potential.  Similarly,  the  problem  of  a  point  charge 
located  between  the  sphere  and  plane,  as  shown  in  figure  3,  will  give 
rise  to  an  infinite  number  of  images.  For  simplicity,  the  series  of 
images  is  broken  into  two  parts  as  shown  in  figure  4a  and  figure  4b . 
The  problem  that  is  treated  here  is  a  charged  conducting  sphere  in 
the  presence  of  an  external  point  charge  and  an  Infinite  ground  plane. 


nn 


•  q' 
-4l 

•  *'l 
Figure  3 
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The  first  series  to  he  treated  is  the  one  resulting  from  the  center 
charge  of  the  sphere  as  shown  in  figure  2.  The  distance  from  sphere  cen¬ 
ter  and  magnitude  of  the  various  charges  can  he  obtained  from  the  simple 
case  considered  in  figure  1.  They  are  given  in  the  following ■ table  to¬ 
gether  with  tjhe  potential  due  to  each  charge  at  the  point  P  located  at 
distance  Z  from  the  sphere  center  as  shown  in  figure  4. 


Magnitude  of  charge 


TABLE  I 

Distance  from  Senter  of  sphere 


Potential  at  P 


The  last  relation  for  can  he  used  to  express  any  change  in  terms  of  Q 
or 

0  =  _2_ 

%  D  D 


n-1 


i  * 


D 


If  we  let  —  =  R  and  —  =  p  then  we  have  R  =  2p 
ana  n 

This  equation  becomes  simplified  if  we  let 


R 


n-l 


with  R^  s=  2p. 


tr  ,  u  ,  u  . 

„  h+l  n+1  .  i  n-l 
Rn  =  T”  °r  ~  ~  U  "  ~ 


n 


n 


9 


then 


(1) 


U  =  2pU  -  U  , 
n+1  *n  n-1 

This  equation  is  a  standard  second  order  differenee^equation.  The 
solution  of  this  equation  can  "be  found  by  assuming  U/  p  where  p  is  to 
be  determined.  If  this  is  substituted  into  equation  1,  we  have 

pn+1  _  2^n  +  pH-1  =  o 

which  has  the  nontrivial  solutions  :  P  =  H  ~  ' -"l.  (2) 

We  then  chose 


Un  =  A1®1  +  k2^2 

where  f3^  =  |i  +  Vp^-l 

*  ,|i.  -  *7^  -  1  j 


n 


U, 


and  the  constants  are  determined  from  the  condition  R^  =  jj~  **  2p. 

u  .  1 
Since  just  the  ratio  ■  P+J~  is  required,  this  one  boundary  condition 

n 

is  sufficient.  The  resulting  solution  is: 

n+1  n+1 

P1  "  P2 
n  =  ‘ 


Pin  -  P2n 


a 

Now  since 


D  D  . 
n  n-1 


—  Q  =  —  — — 

D,  H  R  R  . 
1  n  n-1 


R„ 


Q, 


then 


e,n+1  - 


Q. 


The  total  charge  on  the  sphere  due  to  this  series  of  images  is 


(3) 


—  Q,  +  +  ^2  *  *  * " 

°?  pi  ~ 

-  Q  +  Q  n+1  „  n+1 


1  P1 


-  P, 
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1 


CO 

=  SOj.  -  P2)  I 


and  since  ft.  = 


1  ft- 


o  p1n+1  -  e2n+1  ' 


q(Pl  -  e2)  L 


00  ft 


n 


1  1  -  ft 


2n 


(5) 


is 


The  potential  at  P  due  to  the  series  of  images  inside  the  sphere 


«  oo 

v  =  a  *  y 

%  z  k 


Now 


Z  + 


2 

a 

D 

n 


Q. 


n 


2 

z  +  §- 
n 


Q,  h  ~  P2 
a  z  +  ft. 


L^l 


n+1 


1 

-2TPT 


n+1 


z+Po  2  J 


and 


V 


=  a 

Qrs  a 


P1  -  P2  S 
z  +  P2  i 


ft 


n 


1  !^62n 
"  z+ft-  P2 


(6) 


where  z  =  — 
a 

If  we  set 


n 


@2>  •  i  gD 

i  i-\pQ  n 


then 


Qq  =  Q(ftx  -  ft2)  F(lj^)* 


01 


and 


fti  -  ft,  /  z  +  ft, 


V 


_  q  n. _ _2 


Qfi  a  z  +  P2  "  V  2  +  P2  '  ^2 


(8) 


*  A  Table  of  F(l,ft>)  with  the  first  three  differences  has  been  included 
in  Appendix  A. 
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The  potential  at  P  due  to  the  images  in  the  plane  is 


V, 


% 


s _ r 

+  2h  “ 


oo 


Z  +  2h  -  g- 
n 


and 


Z  +  2h  - 


D 


n 


Q 

a 


-  pg> 

Z  +  P, 


3  n+1  Z  +  ^2  p  n+1  j 
1  '  z  +  p1  P2 


then 


V  =  _  a  ^  F  (  Z  +  B 

%  a  z  +  pi  \ z  +  *V 


(9) 


The  potential  at  P  due  to  all  the  images  is  then 


V 


Q 


=  ^ 


(gl  -  ^2  ) 
Z  +  Po 


Z  +  pr 


z  +  P2 

z  +  p.  '  &2, 


(10) 


A  table  similar  to  Table  I  can  be  made  for  the  series  shown  in 
figure  4a.  See  i?able  II  on  following  page. 

The  same  difference  equation  is  obtained  by  the  substitution 
previously  used^  the  only  difference  being  the  constants  in  the  solu¬ 
tion.  If  the  condition 

R_  =  y,  ( 7-  5),  is  used  then 

.0  a 

Un  '  Pl“  -  Ft!  *2  °r  lettinE  T'  “ 


n 


n 


D  U  _ 
n  n+1 


Un  =  ?i  -  T*  P2  where  again  — • 
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TABLE  II 


Magnitude  of  charge  Distance  from  center  of  sphere  ,  Potential  at  P 


»'o  =  7q 


q'l  =  - g  q\ 

2h  -  & 


q’n  =  r  q’n-l 
n 


2 

2h  -  £- 


D 


n 


4', 


,2 

q  2  "  2  q  1 

2 

2h - 

2h-  a  2 

2h  -  — 

2h  -  £- 

Y 

i  Y 

! 

L 

i 

z  + 


q' 


z  + 


a 


2h  -  -Ty 


q'< 


Z  + 


2h  - 


A 


2 


q’ 


n 


Z  + 


a_ 

D 


n 


where  D  =  2h  -  — 
n  D 


n-1 


The  total  charge  Inside  the  sphere  due  to  these  images  is 


and  as  before 

a  ,  =  -q(i-r) 

^  s 


1 


1-T'P2 


2n 


-q(l-r')  F(Y',p2). 


(11) 
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The  potential  at  P  due  to  this  series  of  charges  can  he  calcu¬ 
lated  by  the  same  method  used  in  deriving  equation  5  and  is 


V 


1  =  where A  h r 


Z+Pj 

Z+Pr 


(12) 


The  potential  at  P  due  to  the  image  charges  in  the  plane  is 


V 


.  =  i  where  B  =  T' 


2+ ft- 


qp  a 


z+p. 


1 


. (13) 


The  total  potential  at  P  due  to  both  series  of  images  is 


V  ,  =  V  ,  +  V  ,  =  a  ' 

q  q'  q  s  a  z+^± 


F(B’pa)  '  Spl  F(A’pa) 


(14) 


The  charge  and  potential  due  to  the  two  series  of  images  shown  in 
figure  4b  can  be  obtained  by  the  same  method  used  in  the  two  cases  given 
above.  The  total  charge  inclosed  in  the  sphere  is: 


a  =  q(i-r)  f(y,P2) 

1  y  ■  P2 

where  Y  =  xs-r-  =  - r—  ■ 

Y'  y  -  P-l 


(15) 


The  potential  due  to  the  image  charges  in  the  sphere  and  the  plane  is 


V 

q 


q  d-r) 

a 


Z+Pp 

F<c<f>2>  -i^F(D’fJ2) 


.  a _ L_ 

a  z+2|j.-y 


(16) 


where 


C 


Ifl 

z+P2^ 


l_z+f52 
■  T '  24-Ej 


and  the  last  term  is  the  poten¬ 


tial  due  to  the  image  of  the  charge,  q,  in  the  plane*  The  total  poten¬ 
tial  at  P  due  to  all  the  charges  and  their  images  is  now  just  the  sum 
of  the  potentials  given  by  equations  l4,  l6;  and  the  potential  due  to 

the  charge  itself,  —  — —  " 

&  7  a  z+y 


The  above  equations  give  the  potential  at  the  point  P  if  some  con¬ 
dition  is  given  to  determine  Q  the  charge  at  the  center  of  the  sphere. 
The  two  cases  we  will  consider  are(l)  when  Q  -  0,  i.e.,  when  the  sphere 
is  grounded,,  and  (ll)  when  the  entire  system  is  neutral,  i.e., 

+  \  +  \  -  -1- 
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Case  I.  Grounded  Sphere 


The  condition  Q  =  0  removes  equation  9  from  the  egression  for 
the  potential  at  P.  This  potential  becomes 

V'  ‘  ■  ViTF  F(B>Po)  -  F(A,fr,>  +  =f- 


z+Pc 


Z+Pr 


i-r 

Z+Pn 


F(D,P0)>  +  J  ~  -  ? 


a  z+y  a  z+2p.-y 


where  D  =  ~  and  C  =  ^  from  equation  (ll),  (12),  and  (l6). 
A  tJ 


Now 


(IT) 


CO 

x  - 

-oo  1-Ap, 


n 

hs  -  - I F  <I^a)  +  ITfc  ’ 


(18) 


then 


v .  aiiT-  V i 

a  z+P  ‘  *— 


°9  P, 


n 


2n  z+3, 


z+fL  oo 

1  X 


n 


1  jrOO  1tB32  2  -00  1-Aft 


2n 


(19) 


Now  when  the  sphere  is  at  a  large  distance  from  the  plane,  32*t:0  so 

that  the  potential  is  given  by  V'  =  SL  P_i - —  1  which  is  the 

h  -*•  oo  a  Ly+z  yz+lj 

potential  at  P  when  there  is  no  plane  present. 

Case  II.  Neutral  System 

The  condition  that  the  entire  system  be  neutral  requires  that 

-q=  q  (l-r)  F(Y,Pg)  +  ft(PrP2)  F(l,32)  -  q(l-Y’)  F(Y',pg). 
This  equation  then  determines  the  charge  at  the  center  of  the  sphere 


Q  =  q 


ti-r )  F(rp2)  -i  -  (i-r)  f(t,p2) 


(Px  -  P2)  f(i,p2) 


(20) 


The  potential  of  the  sphere  is  Equation  (l8  can  be  used  to  write 
equation  (10)  in  a  simpler  form^  hence 
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(a) 


a  -  Pa)  °?  P2n 

VQ  a  a(z  +Ta')' 

1  "  -  .  A 


z+Po  P2 


2n 


The  potential  at  P  becomes 

a  (fi  -  Q  f 

a  (z  + 


v  =  vs  + 


ft 


n 


oo  Z+\  6  2n 

1  '^P2 


a. 


pL:-Hl 

(z+ftx) 


y  PgB  f  ^ 

i*  2n  ^  -,  AO  2n 

~oo  l-Bpg  "°°  I'Ap^ 


wherei  V  =  *  and  Q  is  given  by  equation  (20), 


s  a 


(22) 


The  two  equations  given  (22*  19)  are  good  for  larger  distances  from  the 
plane.  If  the  sphere  comes  close  to  the  plane*  it  is  necessary  to  in¬ 
vert  the  series  as  shown  in  appendix  B  to  obtain  series  which  converge 
rapidly  enough  for  practical  computation .  The  results  are  relatively 
simple  for  the  case  of  the  grounded  sphere  or 


V' 


q  1-V  it 

a  z+pi  .  sa  -n 


»  (  n 

z+e2  2a -Pi  a 


1  e-(^l)»2/q  _  9P 

z+^2  A  cos  0g  -  cosh  ^  (2n+l) 


(19a) 


wh.ere 

it  In  ( -B)  ^  it  In  (-A.)  __J  „  n_  o 

©i  =  — — ' — *-  >  02  =  - ^  and  a  =  In 

5.  RESULTS  AND  DISCUSSION 


In  all  the  results  the  potential  was  calculated  at  the  point  on  the  op¬ 
posite  side  of  the  sphere  at  the  same  distance  from  the  sphere  as  the  charge 
i.e;*  z  =  y.  The  potential  difference  between  the  sphere  and  the  point  was 
computed  i.e.*  V  -  Vg,  and  this  was  normalized  to  the  value  of  V  -  V  at 
infinite  distances  from  the  plane.  In  figure  5  the  normalized  potenlial 
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Figure  5.  Normalized  potential  for  the  case  —  *  1.1 


,  l s 


for  the  case  Y/a  =  1.1  is  plotted  versus  the  distance  from  the  plane  to 
the  sphere's  center,  h/a.  At  great  distances  the  curve  follows  an  inverse 
cubic  hut  as  the  sphere  approaches  the  plane  the  potential  rises  much 
faster.  In  all  the  cases  computed  the  results  are  essentially  the  same 
at  large  distances.  In  figure  6  the  normalized  potential  for  a  number  of 
different  values  of  Y/a  are  plotted  versus  distance  between  the  charge 
and  the  plane.  The  normalized  potential  increases  rapidly  as  Y/a  becomes 
small.  This  is  caused  by  the  small  value  of  V  -  Vs  at  infinite  distances 
from  the  plane  not  by  the  increase  in  V  -  Vs  close  to  the  plane.  A  plot 
of  V  -  Vs  at  infinity  and  at  contact  for  q/a  unity  is  shown  in.  figure  7 
as  a  function  of  Y/a  and  shows  that  the  value  of  V  -  Vs  at  Infinity  de¬ 
creases  more  rapidly  than  when  the  charge  touches  the  plane  for  small 
values  of  Y/a. 
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APPENDIX  A 


TABULATED  VALUES  OF  F(l,x) 


X 

F(l,x) 

A1 

A2 

A3 

0.00 

0.00000  00000 

0.05275  69079 

0.00660  84075 

0.00199  87624 

0.05 

0.05275  69079 

0.05936  53154 

0.00844  71699 

0.00220  35883 

0.10 

0.11212  22233 

0.06781  24853 

0.01065  07582 

0.00273  70204 

0.15 

0.17993  47086 

0.07846  32435 

0.01338  77786 

0.00350  91800 

0.20 

0.25839  79521 

0.09185  10221 

0.01689  69586 

0.00462  88066  . 

0,25 

0.35024  89742 

0.10874  79807 

0.02152  57652 

0.00627  17010 

0.30 

0.45899  69549 

0.13027  37459 

0.02779  74662 

0.00873  23824 

0.35 

0.58927  07008 

0.15807  12121 

0.03652  98486 

0.01252  35835 

0.40 

0.74734  19129 

0.19460  10607 

0.04905  34321 

0.01857  84416 

0.45 

0.94194  29736 

0.24365  44928 

0.06763  18737 

0.02868  85544 

0.50 

1.18559  74664 

0.31128  63665 

;  0.09632  04281 

0.04652  73275 

0.55 

1.49688  38329 

0.40760  67946 

0.14284  77556 

0.08027  03311 

0.60 

1.90449  06275 

0.55045  45502 

0.22311  80867 

0.15012  96943 

0.65 

2.45494  51777 

0.77357  26369 

0.37324  77810 

0.31358  95199 

0.70 

3.22851  78146 

1.14682  04179 

0.68683  73009 

0.76998  88800 

0.75 

4.37533  82325 

1.83365  77188 

1.45682  61809 

2.46195  69250 

0.80 

6.20899  59513 

3.29048  38997 

3.91878  31059 

13.49993J80485 

0.85 

9.49947  98510 

7,20926  70056 

17.41872  11544 

0.90 

|  16.70874  68566 

24.62798  81600 

0.95 

41.33673  50166 

13.39993  70411 

10.11793  29922 

16.11990  36082 

0.96 

54.73667  20577 

23.51787  00333 

26.23783  66004 

88.05181  03388 

0.97 

78.25454  20910 

49.75570  66337 

114.28964  69392 

0.98 

128.01024  87247 

164.04535  35729 

0.99 

292.05560  22976 

0.9995 

8868.79733  5 
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APPENDIX  B 


INVERSION  FORMULAS  FOR  SERIES  EXPANSION 

The  series  used  in  computation  of  the  potential  can  be  converted 
using  complex  variables .  The  method  is  to  multiply  the  series  by  a 
function  having  only  simple  poles  at  the  integers  involved  in  the  sum 
and  with  a  residue  of  l/2rti  at  these  poles. 

Consider 


oo 

1  - 

-oo  1-Xx 


2n 


2  —  ry 

with  x  <  1,  Xx  <1,  and  e  =  x. 


-az 


If  we  multiply  the  quantity 


-2az 


by 


2jtiz 
e  -1 


and  integrate  the 


1-Xe 

resulting  expression  around  a  contour  inclosing  only  the  poles  of 
1 


2itiz 
e  -1 


along  the  real  axis,,  we  obtain  the  given  sum 


-on 

e _ 

.  -2an 
1-Xe 


f  e~az  dz 

J  1-\C-2CB  e2*lz-l 


However,  the  contour  we  shall  choose  to  convert  the  sum  is  shown  in 


This  contour  also  encircles  the  poles  of  the  function  in  the  complex 
plane ,  i.e.,  z  =  ^  i  +  ~~£ct  with  n  <  0 .  Letting 


oo 


p  -  L 


- an 


-2a a 


then 


I 


-oo  -ax 
e 


dx  1  ,1  _ 

„  -7T-r - +  t  p  =  p  +  *a  +  2-rtiR 

-2ax  2rtix  ,  2  K  K  2  o 


-co  1-Xe  -'oo  1-Xe  e  -1 

where  R  is  the  sura  of  the  residues  of  the  function  at  the  poles 


z  =  —  i  +  n  <  0,  and  a  is  the  residue  at  z  =  .  Now  if  we  multiply 

a  2a  o  2a 


-az 


l-\e 


-2az 


by 


l-e 


-2niz 


and  integrate  around  an  appropriate  contour  we  have 


oo 

E 


-On 


-az 


dz 


-2an 


I  ..  ..  -2az  ,  '  -2rtiz 

c-  l-\e  l-e 


-oo  l-\e 

If  the  contour  is  chosen  as  in  figure  B2,  we  have 


f°°  e_0!x  dx  1  _  _  ,  1  „„  . 

J  -2ax  -2*ix  +  2  p  p+  2  o  +  2jaR+ 

-oo  l-\e  l.-re 

where  R+  is  the  sum  of  the  residues  of  the  function  at  z  = 

and  a1  is  the  residue  at  z  =  . 

o  2a 


mti~  ‘lnX 
a  +  2a 


2h 


Now  if  we  add  the  results  of  these  two  contours,  we  have 


1 


co  -ax 


oo  l-\e 


-2ax 


dx 


1-e 


-2rtix 


2itix  . 
e  -1 


+  p  +  -  (aQ  +  a'Q) 


or 


=  2p  +  2«i  [E  +  R  ] 
+ 


/. 


OO  -OCX 

e 


-2ccx 


=  p  +  2rti  [R  +  R  ]  -  i  [a  +  a*  ], 
r  +-£100 


-oo  1-Xe 

Since,  the  integral  on  the  left  is  zero,  we  have 


p  “  |  [a0  +  a'0l  -  2iri  +  R+3 


or 


oo  n 

X 


I 


CO 


-oo  l-\x 


2n 


2a  n/\ 


sin  9  it  r  / 

—  h  (-1)  -  2-7# 


1-cos  0 


a 


'TiC 


sin  9 

+3 


cosh 


a 


-  cos  G 


where  0  =  -  In  X  and  X  >  0. 

a 

The  similar  series  when  X  <  0  can  be  converted  by  the  same  technique 
to  give 


co  n 
x 


*■*  co  -(2n+l)jt^/a  - 

it  +  jt  -  COS  ^ 


I!  ^ 

oo  1-Ax  2a  Jx  ct*T\  1  cos  B  .  cosh  *1  (2n+1) 


To  invert  the  series  for  equation  7  a  different  technique  is  necessary 
from  that  used  in  the  preceding  section.  We  consider 


function  has  a  singularity  like  ^in^x^ 


£  — C-  ■  This 

V  -i  2n 
1  1-x 


near  x  =  1.  We  then  consider 


the  function 
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00 


S  =  £ 


1  1-x 


2n 


In 
In  x 


(1-x)  g 

1 


1  \l-e 
As  before  we  multiply  by 


_  n 
x 


1-x 


2n  2n  In  x. 


apd  integrate  around  a  contour  inclos¬ 


ing  the  real  axis  as  shown  in  figure  B3. 


irtn 

The  poles  along  the  imaginary  axis  are  at:  z  =  -  •- —  • 


2  6 


Then 


,  r o  /  -ax  -ax\ 

is  I  /_e _ e _ \  _ 

2  J__  l  ,  -2ax  2ax  /  2rtix 


dx 


r-oD 


co  \l-e 


-i 


e-aiy  e-aiy 


idy 


1.e-2otly  20iy  /  e-2ty_1 


Jti  . 

-  —  Ro  -  1X1  ^ 


R 


n 


where  R  is  the  residue  of 
n 


-CCS 


-Qfiz' 


1-e 


-2az  2az 


1  ,  nit 

: -  at  z  =-i  — . 

2iuz  ^  Of 


Now  we  multiply  by  - 2xiz  anc^  irrteSrate  around  the  contour  shown  in 


Figure  Bk. 
Then 


1-e 


1  f  °° 

IH 

J  o 


-ax 


-ax 


dx 


.  -2ax  2ax  /  n  -2itix 
vJL-e  /1-e  ■>  oo 


e-ioy  e-ioy 


1-e 


-2iay  2ioy 


idy 

l-e2rty 


-  R  -  Jti  _ 

2  o  T  n 


Sr 


-as 


-az 


where  R  is  the  residue  of 
n 


1-e 


-2az  2as 


1-e 


1  ,  *  nit 

■■■  at  Z  —  1  * 

-2itiz  a 


Adding  the  two  results  give 
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s  = 


oo  -ax 
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■^o  1-e 


-2ax 


dx 


oo 


esc  oy 


o  1-e 


2  ny 


dy  +r  SSSJSi  +  Rj 


ay 


1-e 


2rty  2  o 


r°° 

"  /  2ax 

Jo 


oo 


dx  -  Jti  y  (E  +  R  ) 

i-J  n  n 


Now  R  +  R  =  7— r  and  R  +  R  =0. 

o  o  2«i  n  n 


-L 


oo  -ax 
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oo 
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-2ax 
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cos  oy  1 

Oy  sin  oy 
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The  first  integral  can  be  integrated  to  give 


v. 


oo  /'  -ax  -ax  , 

'  -2ax  "  2ax  )  dx  =  2a  ^ln  2  + 
'a  \l-e 


The  second  integral  can  be  broken  up  into  two  parts  as 


r°°V  1  jA  /_!  cos  qy\ 

Jo  |sin  ay  ay/  yoy  ay  y 


dy 


e2*y-l 


The  last  part  can  be  integrated  to  give 


_1 

2a 


In 


sinh 


a 


1  +  In  2 


The  term  left  can  be  integrated  in  series  to  give 


f"  (  1  .  A  JL  ,  f  (32n-L-l)  ,B  ,3  2„-l 

Jo  \sin  Ciy  ay J  e23t y  ^  2(2n) !  ^  2n-l 


Then  I  /  sinh  — 

S^\Y+21n2+k  (  — 


+  ?  (g211"1-!)  /« 

2 ( 2n )  * 


32n-l) 


2  a 


,2n-l 


where  Y  =  0.5772157  .  (Eulers  const.) 

and  are  Bernoup]_i  numbers.* 


Finally 


l 


n 


1  1-x 


2n 


-  s  {y + taa  - “} ♦ f 


(22n~^-l)  \2  2n-l 

2(2n) !  ^B2n-l'  a 


The  series  involving  the  Bernoulli  numbers  is  an  asymptotic  series  and  can 
only  be  used  when  a  is  very  small,  i.e.,  x  near  1.  (This  is  the  same  re¬ 
sult  obtained  by  Schlomilch  (ref  5)  by  a  different  technique. ) 


*  B.  0.  Pierce,  "A  Short  Table  of  Integrals,"  Ginn  and  Company,  1929 >  3rd 
Ed . ,  p .  90 • 
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REMOVAL  OP  EACH  CARD  WILL  BE  NOTED  ON  INSIDE  BACK  COVER,  AND  REMOVED 
CARDS  WILL  BE  TREATED  AS  REQUIRED  BY  THEIR  SECURITY  CLASSIFICATION- 
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